arXiv:hep-ph/9706347vl 13 Jun 1997 


June 12 1997 
TAUP 2432/97 


PARTON DENSITIES IN A NUCLEON 

A. L. Ayala F £ a ) h )* ** , M. B. Gay Ducati and E. M. Levin 

Institute) de Fisica, Univ. Federal do Rio Grande do Sul 
Caixa Postal 15051, 91501-970 Porto Alegre, RS, BRAZIL 

Instituto de Fisica e Matematica, Univ. Federal de Pelotas 
Campus Universitario, Caixa Postal 354, 96010-900, Pelotas, RS, BRAZIL 

c > HEP Department, School of Physics and Astronomy 

Raymond and Beverly Sadder Faculty of Exact Science 
Tel Aviv University, Tel Aviv 69978, ISRAEL 

T Theory Department, Petersburg Nuclear Physics Institute 
188350, Gatchina, St. Petersburg, RUSSIA 

Abstract: In this paper we re-analyse the situation with the shadowing corrections (SC) in 
QCD for the proton deep inelastic structure functions. We reconsider the Glauber - Mueller 
approach for the SC in deep inelastic scattering (DIS) and suggest a new nonlinear evolution 
equation. We argue that this equation solves the problem of the SC in the wide kinematic 
region where as k = «s xG(x, Q 2 ) < 1. Using the new equation we estimate the 
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1 Introduction. 


The main goal of this paper is to reconsider the whole issue of the shadowing corrections ( 
SC ) to the quark and gluon densities in a nucleon. The motivation to investigate the SC 
for parton distributions comes from some inconsistencies in the interpretation of the present 
HERA data, which we will discuss below. 

The experiment |] shows that the deep inelastic structure function F 2 (x,Q 2 ) has a 
steep behaviour in the small x region (ICC 2 > x > 10~ 5 ), even for very small virtualities. 
Indeed, considering F 2 oc x~ x for small x, the experimental data go from A = 0.15 for 
Q 2 = 0.85 GeV 2 to A = 0.4 for Q 2 = 20 GeV 2 . This steep behaviour is well described 
in perturbative QCD by the DGLAP evolution equations 0-i- The phenomenological 
input, with the quark and gluon distributions at initial virtuality Q 2 = Qq, can be chosen 
at sufficiently low values of Q 2 using the backwards evolution of the experimental data 
in the region of Q 2 sa 4 — 5 GeV 2 . The DGLAP evolution describes the data down to 
Q 2 = 1.5 GeV 2 . 

From the above discussion one can conclude that the parton cascade is a rather deluted 
system of partons with small parton - parton interaction, which can be neglected in a first 
approximation. Therefore, no SC is needed to describe the experimental data. 

Nevertheless, there is a set of several facts that does not fit into this scheme. 

1. We expect that the quark distribution will not grow indefinitely as x goes to zero, 
since it would violate unitarity for some value of x ||. In Ref. [|j] the unitarity bound (UB) 
for F 2 (x, Q 2 ) was established and it turns out that the F 2 structure function reaches the UB 
at Q 2 = 1 — 2 GeV 2 in HERA kinematic region. It means that large shadowing corrections 
to the normal DGLAP evolution should take place in this kinematic region. 

2. From HERA data we can evaluate also the probability k of the parton - parton (gluon 
- gluon) interaction, which is given by ||, [^] 


K 


xG(x, Q 2 ) 


a(GG) 

Q 2 nR 2 


3 7r dg 
Q 2 R 2 


xG(x, Q 2 ) , 


( 1 ) 


where xG(x, Q 2 ) is the number of partons ( gluons) in the parton cascade and R 2 is the radius 
of the area populated by gluons in a nucleon. a(GG) is the gluon cross section inside the 
parton cascade and was evaluated in 0 . Using HERA data on photoproduction of J/T meson 
[|I] the value of R 2 was estimated as R 2 < 5 GeV~ 2 ||. Using the GRV parameterization 
[ffj for the gluon structure function and the value of R 2 = 5 GeU~ 2 , we obtain that n reaches 
1 at HERA kinematic region ( see Fig-Ill )> meaning shadowing corrections should not be 
neglected (see Ref. [|5] for details). 
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Figure 1: Contour plot for k = cte for the GRV 95\ f|/ distribution and R 2 = 5 GeV 2 . The 
solid lines show the HERA kinematic region 43/- 

3. The situation looks even more controversial if we plot the average value of the anoma¬ 
lous dimension, < 7 > = d\n[xG(x, Q 2 ))/d\nQ 2 , in the GRV parameterization. Fig.|2l 
shows two remarkable lines: the line < 7 >= 1 , where the deep inelastic cross section reaches 
the value compatible with the geometrical size of the proton; and the line < 7 >= 1 / 2 , 
which is the characteristic line in which vicinity both the BFKL Pomeron [Q and the GLR 
equation || should take over the DGLAP evolution equations. HERA data passed over the 
second line and even for sufficiently small values of Q 2 they crossed the first one without any 
indication of a strange behaviour near these lines. 

From the above discussion we conclude that the physical interpretation of HERA data 
looks controversial and the statement that the DGLAP evolution works is first but not the 
last outcome of HERA data. We expect also that the SC will be essential in HERA kinematic 
region. However, to treat the SC we need a more general approach than the GLR one 0 . 0 , 
which gave a theoretical approach to the SC only to the right of the critical line < 7 > = |, 
at large values of Q 2 . In this paper we develop a general approach to the SC, that will 
allow us to give reliable estimates for the SC in the kinematic region to the left of the line 
< 7 >= 1 / 2 , i.e., for large ln( 1/x) and small Q 2 . 

This paper is organized as follows. In section 2 we present the Glauber or Eikonal 
approach for high energy scattering off nucleon. We discuss the SC for the gluon distribution 
in a nucleon and for the F 2 deep inelastic structure function, concentrating mostly on the 
theory status of this approach and on the predicted scale for the SC. We analyse also the 
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Figure 2: Contours for <7 > = 1 and < 7 > =1/2 for the GRV95 gluon density and 
HERA kinematic region. 

effect of the SC on the ^-dependence of the gluon distributions and F 2 , and the consequent 
approach to the “soft” physics. After a brief discussion, in section 3, of the first correction 
to the Glauber (Eikonal) approach, we present, in section 4, a new evolution equation that 
sums all SC diagrams which are of order (as ln(l/x) In Q 2 ) n and K n . This equation is solved 
for as constant in the semiclassical approach. The solution is compared with that of the 
GLR equation || and with the Glauber approach. Finally, section 5 summarizes our results. 


2 The Eikonal approach in QCD . 

2.1 The Mueller formula. 


In this section we present the general idea for the Glauber formula in QCD both for the 
gluon distribution and the F 2 structure function of DIS. These Glauber formulas will take 
into account the SC from the eikonalization of the parton-nucleon cross section. The ideas 
were originally formulated in Refs.0]- [ITT . 


Let us consider a high energy virtual particle that probes the nucleon gluon distribution 
( for example, such a probe could be a virtual graviton or a heavy Higgs boson). In the 
space-time picture of the process, the virtual probe decays in a gluon-gluon (GG) pair, with 
transverse separation r± and the fraction of energy of the probe z and 1 — z, which interacts 
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with the nucleon through a gluon ladder exchange. In high energy (small x) limit, r± is 
constant during the interaction for £ -C 1/(2 mR), where R is the target nucleon size. 
Therefore, the transverse distance r± is a good degree of freedom for our problem. 

The cross section of the absorption of gluon(G*) with virtuality Q 2 and Bjorken x can 
be written in the form 

r i r ( R r . r r /2 r h 

cr(G*) = / dz / - / --I v b < ?*(Q 2 , r±, x, z) I 2 <7 to t(GG + nucleon ) , (2) 

JO J 7T J 7T 

where z is the fraction of energy carried by the gluon, b t is the impact parameter and 
is the wave function of the transverse polarized gluon in the virtual probe. The cross 
section of the interaction of the GG pair with the nucleon is a to t(GG + nucleon ) = a(x, r\). 
It is important to note that this description is valid in leading log ( 1/x ) approximation 
(LL(l/x)A), in which only contributions of the order of ( ctg ln(- )) n were taken into account. 
In this approximation we can neglect the change of z during the interaction and describe the 
cross section a to t(GG + nucleon ) as a function of the variable x 0, 

The cross section cr(a:, 7j_) can be written in the form 

a(x,r 2 L ) =2 J d 2 b t Ima(x,r±,b t ) , (3) 

where a is the elastic amplitude for which we have the s - channel unitarity constraint 

2 Ima(x, r j_, b t ) = \a(x,r ± ,b t )\ 2 + G in (x, r x , b t ) , (4) 

where G in is the contribution of all inelastic processes. Eq. (J3j) gives the relation between two 
unknowns a and Gi n and has a general solution in the limit of small real part of the elastic 
amplitude at x —» 0 


m, [i3 


a(x, r±,b t ) = i { 1 — e h n G,rxM) j. ^ ( 5 ) 

G in (x,r ± ,b t ) = 1 - 

The opacity function Q has a simple physical meaning, namely e ' is the probability 
that a GG pair has no inelastic interaction during the passage through the target, is an 
arbitrary real function, which can be specified only in a more detailed theory or approach 
than the unitarity constraint. One of such specific model is the Glauber approach or the 
Eikonal model. 

However, before we discuss this model let us make one important remark on the strategy 
for the approach to the SC: being able to calculate the opacity 0 we will have the theory or 
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model for all inelastic processes. Indeed, using AGK cutting rules |I5| we can calculate any 
inelastic process, if we know in accordance with the s-channel unitarity. It is worthwhile 
mentioning that the inverse procedure does not work. If we know the SC in all details for a 
particular inelastic process, say for the inclusive production, we cannot reconstruct all other 
process and the total cross section in particular. 

Now, let us built the Glauber approach. First, let us assume that is small (D -C 
1) and its b t dependence can be factorized as Q = Q(x,r±) S(b t ) with the normalization 
/ d 2 btS(b t ) = 1. The factorization was proven for the DGLAP evolution equations § 
and, therefore, all our further calculations will be valid for the DGLAP evolution equations 
|I1 in the region of small x or, in other words, in the Double Log Approximation (DLA) 
of perturbative QCD (pQCD). Expanding Eq. (^]) and substituting it in Eq. (2|), one can 
obtain 


<Jtot(x : r_ l) = D(x,r ± ). 

The substitution of Eq. (g) in Eq. (g) gives 


a(G*) = 


dz 


d 2 r\ 


7T 


Vf(Q 2 ,r ± ,x,z ) | 2 f2(x, r_|_) 


( 6 ) 


(7) 


Using Eq. (|7|) and cx(G*) = xG(x , Q 2 ) as well as the expression for the wave function 


of the GG pair in the virtual gluon probe 
the glnon structure function P-flOl 


see Ref. flT0| ), one can express D(a:, r±) through 


n = 


a 


GG, 
N 1 


X, 


3 7r 2 as 


r 2 xG(x, — 2 ~) 


( 8 ) 


The Glauber (eikonal ) approach is the assumption that D = i}S(b t ) with D given by 
Eq. (HD in the whole kinematical region. From the point of view of the structure of the final 
state, this assumption means that the typical rich inelastic event was modeled as a sum of 
the diffractive dissociation of the GG pair and uniform in rapidity distribution of produced 
gluons. For example, we neglected in the Glauber approach all the rich structure of the large 
rapidity gaps events including the diffractive dissociation in the region of large mass. 

Substituting Eq. (|8|) in Eq. (|) and the result in Eq. (Q), and using the wave function 
r _|_, x, z) calculated in Ref.[jl(| we obtain the Glauber (Mueller) formula for the gluon 
structure function 


xG(x , Q 2 ) 



<7%° (x' 


(9) 
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The first term in the expansion of Eq. (J9|) with respect to cr GG gives the DGLAP equation 
in the region of small x. 

To simplify our calculations we will use the Gaussian parameterization for the profile 
function S(b t ), namely 

1 


S(bf) = 


7 TR 2 


e r z 


( 10 ) 


Using this profile function and integrating over b t , we obtain (N c = Nf = 3) 
2 R 2 dx' h dr 2 . 


xG(x, Q 2 ) = — y - f - 4 - / Q5 { C + ln(n G (x\rl)) + Ei{k g (x', r ^))} , ( 11 ) 

7T J x X J a T“ | 


where C is the Euler constant and £1 is the exponential integral (see Ref.flql Eq. 5 . 1 . 11 ) 
and 

3as(4-)7r?"^ 1 

M*', r|) = x > G dglap { x , (12) 


The Eq. ([|) is the master equation of this section and it gives a way to estimate the value 
of the SC. 

A similar Glauber (eikonal ) formula may be obtained for the deep inelastic structure 
function F 2 (x, Q 2 )- In this case, the virtual probe decays into a quark anti-quark pair which 
interacts with the nucleon trough a gluon ladder exchange. Taking into account Nf quark 
flavours and integrating the quark pair wave function over z, we obtain a Glauber formula 
for F 2 ( see Refs. [|] || |K| ) 


N r . Nf 


f 2 (x,q 2 ) = -y z z 

0 IT 6 i 


r°° dr 2 . 


J d 2 b t { 1 


g- \ ^ qq(x,r±,b t ) 


}. 


(13) 


where S{b 2 ) = |h2, and Zf is the quark charge fraction. After integration over 

b t using the Gaussian profile function S(b t ) one can reduce Eq. ([Hip to ( N c = 3 ) 


JD 2 Nf i- •*■ ,y 

f 2 (x,q 2 ) = 7^J2z 2 f J-pr{ c + + E i K^'ui))} i ( 14 ) 


-ki dr 


27r 2 | 


with = 4 kq- 


2.2 Theory status of the Mueller formula. 

In this section we shall recall the main assumptions that have been made to obtain the 
Mueller formula. 
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1. The gluon energy (x) should be high (small) enough to satisfy x -C ^ and 
otsln(l/x) < 1. The last condition means that we have to assume the leading ln(l/x) 
approximation of perturbative QCD for the nucleon gluon structure function. 

2. The DGLAP evolution equations hold in the region of small x or, in other words, 
a.$ln{l/r \) < 1. One of the lessons from HERA data is the fact that the DGLAP evolution 
can describe the experimental data. 


These two assumptions mean that we describe the gluon emission in Double Log Ap¬ 
proximation ( DLA) of perturbative QCD. In other words, we extract from each Feynman 
diagram of order ctg the contribution of order (as lnl/x lnQ 2 /Q(j ) n , neglecting all other 
contributions of the same diagram. In terms of the DGLAP evolution, we have to assume 
that the DGLAP evolution equations describe the gluon emission in the region of small x. 
However, the first assumption is very important for the whole picture, since it allows us to 
treat successive rescatterings as independent and simplifies all formulae reducing the problem 
to an eikonal picture of the classical propagation of a relativistic particle with high energy 
(ER 1) through the target. 

3. Only the fastest partons (GG pairs) interact with the target. This assumption is an 
artifact of the Glauber approach, which looks strange in the parton picture of the interaction. 
Indeed, in the parton model we rather expect that all partons, not only the fastest one, should 
interact with the target. In the next section we will show that corrections to the Glauber 
approach due to the interaction of slower partons are essential in QCD. 


4. There are no correlations (interaction) between partons from the different parton 
cascades. This assumption means that even the interaction of the fastest GG-pair was taken 
into account in the Mueller formula only approximately and we have to assume that we 
are dealing with large number of colours to trust the Mueller formula. Indeed, it has been 
proven that correlations between partons from different parton cascades lead to corrections 
to the Mueller formula of the order of 1/iV 2 , where N c is the number of colours (||18||). 


Let us discuss the large N c approximation in more detail. The main principle, used in our 
large N c approach, has been formulated by Veneziano et al. 0 , namely, we sum leading N c 
diagrams considering N c a§ ~ 1 while ag < 1, separately in each topological configuration. 
For example, at high energy ( low x ) the main contribution gives the planar diagrams ( see 
Fig#. )which lead to the DGLAP evolution equations in leading ln(l/x) approximation of 
pQCD. In fact, it has not been proven yet that more complicated planar diagrams such as of 
Fig# could be neglected at high energies, but several simple examples have been considered 
( see Ref.Jl8j ) which show that this hypothesis is very likely ( see review [BJ for a detailed 
discussion on the subject). The two sheet configuration (see Fig#) is suppressed since it 
is proportional to an extra as- The main contribution at high energy in this configuration 
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are the diagrams that has been summed in the Glauber - Mueller approach (see Ref. ||18||) in 
which the smallness of the order of cts is compensated by extra power of xG(x, Q 2 ). The 
actual parameter of the Glauber - Mueller approach is k which could be of the order of 1. 
However we neglected the contribution of the diagrams of Fig. 0 d type in which one can see 
a transition from one two sheet topology to another. As has been shown |n| such transition 
is suppressed by a factor -Tj. The use of the A- expansion can be suspicious in this case, 
because such transition gives a corrections to the anomalous dimension. However, it was 
shown |25| that A- - expansion works quite well, justifying our approach. To be safe we 
assume that ln(l/x) < 1 . 


5. There are no correlations between different gluons inside the nucleon, of about the 
same energy, except the fact that they are distributed in the area with radius R. This is 
the main assumption of the Eikonal model which makes this model simple and, of course, 
it is a great simplification of the unknown ( nonperturbative ) structure of the nucleon. R 
is the correlation radius for such gluons and the estimates from the HERA data have been 
discussed in the introduction. 


2.3 The modified Mueller formula. 


The next step of our approach is to give an estimate of the SC using the Mueller formula. 
However, before doing so, we have to study how well the DLA of perturbative QCD works, 
as it was heavily used in the derivation of the Mueller formula. Let us consider the small x 
limit of the DGLAP evolution equation, which takes a simple form 


d 2 xG(x, Q 2 ) «s 
d\n(l/x) d\nQ 2 n 


xG(x, Q 2 ) . 


(15) 


This expression may be obtained from the full DGLAP equation taking only the DLA limit 
of the anomalous dimension, namely, 7 DLA = as j Vc -. 

In order to estimate numericaly how well the DLA works, we use the GRV parameter¬ 
ization |Q] for the nucleon gluon distribution. This parameterization describes all available 
experimental data quite well, including recent HERA data at low x and low Q 2 ( Q 2 > 
1.5 GeV 2 ). Moreover, the GRV parameterization is suited for our purpose because (i) the 
initial virtuality for the DGLAP evolution is small (Q q « 0.25 GeV 2 ) and we can discuss the 
contribution of the large distances in the MF having some support from experimental data; 
(ii) in this parameterization the most essential contribution comes from the region where 
a s lnQ 2 ~ 1 and a s lnl/x ~ 1. This allows the use of the double leading log approximation 


of pQCD, where the MF is proven [10]. It should be also stressed that we consider the GRV 
























Figure 3: Topological configurations for large N c approach: a) ladder diagram; b) general 
planar diagram; c) Glauber-Mueller approach; d) general two sheet diagram. 
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parameterization as a solution of the DGLAP evolution equations, disregarding how much 
of the SC has been taken into account in this parameterization in the form of the initial 
gluon distribution. 


However, in spite of the fact that the DGLAP evolution in the GRV parameterization 
starts from very low virtuality ( Qq ~ 0.25 GeV 2 ) it does not reach the DLA in the accessible 
kinematic region ( Q 2 > 1 GeV 2 , x > 10~ 5 ). To illustrate this statement we plot in Fig .[3] 
the ratio: 


< a s > 


GRV 


a s 


d 2 xG GRV (x,Q 2 ) 


d\n(\/x) 9 In Q 2 

^xG GRV (x ,Q2) ' 


This ratio is equal to 1 if the DLA holds. From Fig.|] one can see that this ratio is rather 
around 1/2 even at large values of Q 2 . 



Figure 4: The ration <ag / > s 


GRV 


for different values of Q 2 in the GRV95 parameterization. 


We can understand why the corrections to the DLA are so big modeling the small x 


anomalous dimension 7 by a simple formula [26 


7M = ^{I- 1} , (16) 

7r u 

which has the correct DLA limit at small u and satisfies the momentum conservation con¬ 
straint 7 ( 0 ; = 1) = 0. The typical values of uj in all available parameterizations, even in the 
GRV, which is the closest one to the DLA, is < u > ~ 0.5. Therefore, there is about 50% 
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correction to the DLA, meaning it cannot provide reliable estimates for the gluon structure 
function. 


On the other hand, our master equation (see Eq. (|9])) is proven in DLA. Willing to 
develop a realistic approach in the region of not ultra small x {x > ICE 4 ) we have to change 
the master equation ( Eq. (|9|) ). We suggest to substitute the full DGLAP kernel in the first 
term of the r.li.s, which gives 


xG(x, Q 2 ) 


xG(x,Q 2 )(Eq. ( 9)) + xG dglap (x,Q 2 
a s N c r 1 rQ 2 dx' dQ' 2 

E JQl 


— ^x'G 


DGLAP, 


7T 


X' 


Q 


12 


x', Q 


/ 2 \ 


(17) 


The above equation includes xG DGLAP (x , Qq) as the initial condition for the gluon distri¬ 
bution and gives xG DGLAP (x,Q 2 ) as the first term of the expansion with respect to kg- 
Therefore, this equation is an attempt to include the full expression for the anomalous di¬ 
mension for the scattering off each nucleon, while we use the DLA to take into account all 
SC. Our hope, which we will confirm by numerical calculation, is that the SC are small 
enough for x > 10” 3 and we can be not so careful in the accuracy of their calculation in 
this kinematic region. Going to smaller x, the DLA becomes better and Eq. © tends to 
the master equation 

We suggest an analogous improvement in the calculation of the SC correction for the 
structure function F 2 . We subtract the Born term of Eq. ( |T4|) and add the DGLAP evolved 
structure function F PGLAP . The result reads 


F2(x y Q 2 ) 


pSC _ pBORN pDGLAP 


(18) 


where F pc is given by Eq. (jL|), and the Born term is given by 


pBORN ( 


Nf 

E Q 2 ) = 2 ^ 12 «s z 


/•In Q 2 2 
An Q g 3 


d{lnQ u )xG 


-i DGLAP, 


X, Q 


/ 2 \ 


(19) 


As discussed above, F PGLAP (x , Qq) is the initial condition for F 2 and F PGLAP (x , Q 2 ) is the 
hrst term of the expansion with respect to K q . Now, we are able to evaluate the amount 
of SC in the nucleon gluon distribution and in the F 2 structure function predicted by the 
eikonal approach. 


2.4 The gluon structure function for nucleon. 

Before performing calculations we would like to discuss which observable we are able to 
calculate using the Mueller formula. One can see from Eq.(9) that the large distances enter 
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into this formula. In spite of the fact that this formula provides the infrared safe answer, 
or in other words, the integral is finite in the region of the large distances, we cannot trust 
our calculations in this region since we have no arguments to apply pQCD in this kinematic 
region. Therefore, first we need to indroduce a scale of distances for which we can trust the 
pQCD calculations. We decided to choose the value of r\ 0 = 1 GeV~ 2 as such scale. The 
argument comes from the new HERA data on the deep inelastic structure function F 2 (^', Q 2 ) 
[I] which show that the GRV parameterization works quite well starting from Q 2 > 1 GeV 2 . 
Therefore, we consider the GRV structure functions at Q 2 = 1 GeV 2 both for gluons and 
quarks, as a lucky guess of the initial conditions for the QCD evolution that takes into account 
the unknown nonperturbative large distance physics, disregarding all physical motivation for 
such a guess. We believe that such a choice give us a reliable estimate of the SC which stem 
from small distances, being on the theoretical control of pQCD. In the following calculations 
we will use a fixed value for the QCD running coupling constant, cts = 0.25. 

In order to evaluate quantitatively the influence of the SC on the evolution of gluon 
structure function we calculate three functions 


xG(x, Q 2 ) 
xG GRV (x,Q 2 ) ’ 

9 In (xG(x, Q 2 )) 
d\n(l/x) 

d\n(xG(x, Q 2 )) 
1 01n(Q7Q§) 


( 20 ) 

( 21 ) 

( 22 ) 


We chose these functions because, in the semiclassical approach ( see Ref. || ), the gluon 
structure function has a general form 


xG(x,Q 2 ) = C{Q 2 } <7> {-} <U}> 

x 


(23) 


where we expect the exponents to be smooth functions of In(-) and ln(Q 2 ). We will see 
a posteriori that, indeed, < to > and < 7 > turn out to be smooth functions in the 
region of small x. All functions in Eq. (|20|) - Eq. (22) have transparent physical mean¬ 
ing. Indeed, the function < 7 > is the average anomalous dimension in DIS, the func¬ 
tion < u > is the intercept of so called “hard” Pomeron while function R\ is equal to 
l/Np orneron , where Np omeron is the average number of “hard” Pomerons taking part in the 
rescatterings . “Hard” Pomeron in our approach is xG(x, -r) in Eq. @ ( (xG(x, 4-)) n cor¬ 
responds to n - Pomeron exchange). Indeed, due to the AGK cutting rules [Hj, which 
have been proven in QCD in Ref. [E7|, the inclusive cross section is proportional to the 
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Figure 5: R\ ratio as a function of y = ln(l/x) and Q 2 for two values of R 2 . 

exchange of one Pomeron or, in other words, only the first term in Eq. (jD[) contributes to 
the inclusive cross section of particles in the central region of rapidity (ai nc i = ^). On 
the other hand, a inc i = N a tot , where N is the average multiplicity and a tot is the total 
cross section given by Eq. H). At is equal to N Pomeron N ( one Pomeron exchange ), where 

Nione Pomeron exchange ) = --- filTp ° mer °" . Noticing that <7p omeron is the first term in 

Eq. we obtain that N Pomeron = 

In Fig.|5|, we show the values of R\ as a function of y = ln{l/x) for different values of 
the virtuality Q 2 and two values of R 2 . The hgure illustrates the general feature of the SC 
in the eikonal model: the distribution xG is strongly suppressed as y increases ( x tends to 
zero). In HERA kinematic region, 2 < y < 10 and Q 2 = 2.5 GeV 2 , the modified Mueller 
formula for R 2 = 10 GeV~ 2 predicts a suppression that varies from less then 1% up to 10%. 
For Q 2 = 20 GeV 2 , the suppression varies from less then 1% up to 13%. This effect is even 
bigger for smaller value of R 2 . For example, for R 2 = 5 GeV~ 2 and Q 2 = 2.5 GeV 2 , the 
suppression goes from less then 1% for y = 2 to almost 15% for y = 10. For Q 2 = 20 GeV 2 , 
the suppression goes from less then 1% for y — 2 to 22% for y = 10. 

We can interpretate the result shown in Fig.|| in a different way saying that the average 
number of the “hard” Pomerons participating in rescattering changes from 1 at x = 10 -2 to 
1.25 (R 2 = 10 GeV~ 2 and Q 2 = 10 GeV 2 ) or 1.43 (R 2 = 5 GeV~ 2 ) at x ^ 10" 5 . 

In Fig.|] are shown the results for < u > as a function of y for several values of Q 2 and 
two values of R 2 . As we can see, the values of u are softened by the SC when compared 
with the GRV values. 
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Figure 6: The values of < u > as a function of y = ln( 1/x) and Q 2 for GRV and two values 
of K 2 . It is also plotted the soft Pomeron intercept u> = 0.08. 
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In Fig.|7| are shown the results for < 7 > calculated from expression (p2|). The values 
of < 7 > suffer a strong suppression as compared to the GRV ones. It means that the Q 2 
evolution of the gluon distribution is significantly modified by the SC. 

The figures § | and 0 also illustrate the dependence of SC on the value of R 2 . For 
example, for R 2 = 10 GeV~ 2 , the SC, generated by the eikonal approach, suppress the 
value of ui, but the value of the SC is not enough to provide a match between the “soft” 
and “hard” processes. To illustrate this point we put in Fig.| 6 ] the intercept of the so called 


28 . For R 2 


“soft” Pomeron a i so f t = 0.08 
does not reach the “soft” Pomeron intercept 


= 5 GeV~ 2 , the suppression is bigger, but u 
However, for Q 2 = 20 GeV 2 sufficiently small 
values of x (y ~ 20) and Q 2 = 20 GeV 2 , the SC reduces the values of u reaching the “soft” 
Pomeron intercept. One can see that < a; > turns out to be bigger than this intercept at 
small values of Q 2 . It occurs because these values of Q 2 are close to the initial condition and 
do not generate enough shadowing. 


As we have discussed the SC in the eikonal approach considerably reduce the value of 
< 7 >. For R 2 = 5 GeV~ 2 , the mean value of the anomalous dimension < 7 > reaches the 
value | only for Q 2 = 2.5 GeV 2 in the range 10 < y < 16. For Q 2 > 2.5 GeV 2 , < 7 > is 
always less then |. Thus we expect that the BFKL Pomeron will not be seen in the deep 
inelastic gluon structure function at least in HERA kinematic region. For R 2 = 10 GeV~ 2 , 
however, the strength of the SC is such that < 7 > > in the accessible kinematic region. 
It means that, for Q 2 < 2.5 GeV 2 and 7.5 < y, the BFKL Pomeron could give a visible 
contribution. The above discussion illustrates the role of SC and the R 2 parameter on the 
transition to the BFKL dynamics which we are going to discuss a bit later in more details. 

In order to complete the discussion of the SC from Mueller formula, we plot in Fig.|8| the 
expression 


< off > 

< a s > GRV 


xG GRV {x,Q 2 ) 

9 (tG -It r> 2 A 


dln±dlnQ 2 V )) 

X ^ 


xG n (x, Q 2 ) 


9 (rr-nCRVi^ niw 

ainiSln ^ )> 

X ^ J 


(24) 


which is the ratio of the effective as in DLA with and without the SC. In some sense this 
is the effective value of the QCD coupling constant in the parton cascade with SC, and 
characterizes the amount of shadowing taken into account in the DGLAP evolution of the 
gluon distribution. The results present a weakening of the effective coupling constant which 
accepts a simple interpretation. In the small-x region, the gluon density increases and in a 
typical parton cascade cell there are many partons with different colors. The color average 
gives smaller mean color charge at large gluon densities and < as > is reduced. From the 
figure we can see that SC modify the DGLAP evolution for a large kinematic region, even 
for Q 2 > 10 GeV 2 , though including a substantial part of HERA data. 
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Figure 8 : The values of the effective 
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«s in the parton cascade with SC. 


2.5 The DIS structure function F 2 . 


In this section we evaluate the SC predicted by the Glauber approach for F 2 . We calculate 
F ‘2 from the modified Mueller formula (IB) 


rpBORN 1 t^DGLAP 
r 2 ' r 2 


F 2 (x,Q 2 ) = F 2 b ‘ 

where F 2 C and F 2 ORN are given by Eq. O and Eq. flT9|) respectively. We calculate F 2 GLAP 
from the GRV 95 |29j] distribution u, d, s and xG using the expression 


F. 


DGLAP 


(x, Q 2 ) = J2 4 x [<l( x > Q 2 ) + q{x, Q 2 )] + T 2 c O, Q 2 , m 2 c ), 

u,d,s 


(25) 


where q and q represent the light quarks distributions. F 2 describes the heavy flavour 
contribution to F 2 (disregarding bottom contribution). I 11 this scheme the heavy quarks 
are not considered as intrinsic partons in the Renormalization Group evolution, but are 
produced, in LO of perturbation theory, by the gluon- 7 * fusion process. Thus, the charm 
component F 2 will be generated perturbatively from the intrinsic gluon distribution. 


I 11 Fig.|9| we plot F 2 using the full DGLAP evolution, given by Eq. (|25D , and the modified 
Mueller formula (f 

the gluon distribution in LO, as described in Ref. 


to take into account the SC. The F 2 component was calculated from 


distribution that enter in the F 2 calculation. We expect that it would give a correction for 


F 2 smaller then the correction given by Eq. 


9|. The SC are not included in the gluon 

a 

, at least by a factor as(4m^) 
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Figure 9: F 2 calculated from fl 7$) for Q 2 = 1.0 GeV 2 . 

In figs. ®) and (|ll|) we present the calculation of Ff )CLAP an d F 2 given by Eq. ( ]18|) as 
a function of x for different values of Q 2 . The results are compared with HI 96 data. As 
we can see from the figures, the SC are small for Q 2 = 1.5 GeV 2 since it is very close to 
the initial value Ql = 1 GeV 2 . When Q 2 grows from Q 2 = 1.5 GeV 2 to Q 2 = 6 GeV 2 , the 
SC increase. For Q 2 > 6 GeV 2 , the SC start to fall down again. In fact, the figures show 
the amount of shadowing that we would expect to be present in the experimental data, if 
the Glauber approach is supposed to be correct. We conclude that the contribution of the 
SC to F 2 (x,Q 2 ), predicted from the eikonal approach are rather small in HERA kinematic 
region unlike the case of the gluon structure function that has been discussed in the previous 
section. This is an explanation why we are able to describe the HERA experimental data 
on F 2 (x,Q 2 ) in the GLAP evolution, without taking into account SC in spite of the fact 
that parameter n is rather big. It turns out the story of the SC contribution is the story 
about the gluon densities for which we have no independent experimental information. The 
value of Ri for xG(x, Q 2 ) is rather big, as we have discussed, but it is still smaller than 
uncertainties in the values of xG(x, Q 2 ) from different attempts ( GRV, MRS and CTEQ) 
to describe the experimental HERA data in the framework of the DGLAP evolution, that 
exist on the market. 
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Figure 10: F 2 calculated from fll\ ) for Q q = 1.0 GeV 2 compared with HI96 dataf^j. 
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Figure 11: F 2 calculated from ([ T$) for Ql = 1.0 GeV 2 compared with H196 datafl ]/. 
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2.6 The b t dependence of the gluon distribution. 


In the next two subsections we will investigate how the SC work at different values of the 
impact parameter b t . We hope, that the impact parameter distribution will help us to 
understand better the matching between “soft” and “hard” interactions. In order to do this 
we introduce the function xG(x, Q 2 , b t ) defined by the normalization 

J d 2 b t xG(x, Q 2 , b t ) = xG(x, Q 2 ) , (26) 

where xG(x, Q 2 ) is the gluon distribution in perturbative QCD. We introduce also the func¬ 
tion xG(x, Q 2 , t ) given by 

xG(x,Q 2 ,t) = J d 2 b t e l ^ t xG(x,Q 2 ,b t ) . (27) 


where t is the momentum transfer t = —q 2 . 

In the DGLAP evolution equation, the b t dependence can be factorized out in the form 

0 

xG{x,Q 2 M ) = xG(x,Q 2 )S(b t ) , (28) 

where the profile function S(bt) is related to the two gluon nucleon form factor F(t) as 


F(t) = J dV*® * 1 S(b t ) . 


(29) 


If we consider the exponential parameterization of the profile function S(bt ), Eq.( ^) and 
( |2U| ) give also an exponential t- dependence for the gluon structure function, which can be 
obtained by the ratio 


xG(x, Q 2 , t) 
xG(x, Q 2 , t — 0) 



(30) 


Thus the DGLAP evolution predicts a slope B 


/P 

4 ’ 


that does not depend on x and Q 2 . 


However, as we will show, the SC will lead to an x and Q 2 - dependence of the B slope, 
as well as a non factorizable b t - dependence of the gluon structure function. 


We introduce the impact parameter b t dependent function xG(x, Q 2 , b t ) defined by Eq. (|26|), 
where xG(x,Q 2 ) is the gluon distribution given by the Glauber formula ([]). From (|?6|), 
xG(x, Q 2 , b t ) can be written as 


~ 2 , v f 1 dx’ [°° dr\ 

xG(x, Q ,b t ) = — — — j- { 1 

7C Jx x ' JAr r. 


i „gg 


S(bt) 


}, 


(31) 


21 





where a GG was defined in Eq. (| 8 [). Following the discussion before Eq. ([L7|) we will modify 
Eq. ([H]) in order to use the full DGLAP kernel to describe the gluon-nucleon interaction and 
the DLA to take into account the SC. We add a term proportional to the gluon distribution 


and subtract the Born term of ( fTT] ) . This procedure gives 


xG(x, Q 2 , b t ) = 


xG 


7T 


DGLAP , 


x,Q 2 )S(b t ) + xG(x,Q 2 ,b t )( 31) 


r l dx' 


dr] 


x' 


-rr GG 
2 N 


S(b t ). 


(32) 


This expression keeps the same b t dependence in the Born term and gives the correct match 
with Eq. (26). In all calculations we use a fixed value for the coupling constant «s = 0.25 
and put Ql = 1 GeV 2 . 


In order to investigate the general properties of the b t dependent distribution, we calculate 
all the three functions Ri (Eq. (|20|) ), < u > (Eq. (|2iP ) and < 7 > (Eq. (|22|) ) using the above 
expression for xG(x,Q 2 ,b t ) instead of xG(x,Q 2 ). The results are plotted in Fig. [12] for 
b t = 0. We can see from the figure that the SC are bigger for the b 2 = 0 case when 
compared to the b ? integrated results. The correction term in (|32|) will be proportional to 


{1 — kg — gxp(-kg)}, which is bigger then the nonlinear term of Eq. (|l7j). 


It is interesting to notice that at b t = 0 the SC give < u > approaching a j so f t = 0.08 for 
y ps 14, which corresponds to x = 10 - '. Thus we may expect eventually a matching of the 
DIS data with the “soft” processes at this very small value of x. 

The b t dependence of the gluon distribution is illustrated also in Fig.[L3|, which shows 
xG(x,Q 2 ,b t ) normalized to b t = 0 at different values of x, Q 2 and R 2 = 5 GeV~ 2 . Due to 
the SC, the width of the b t - distribution becomes bigger with respect to the exponential 
profile function, and presents an x and Q 2 dependence. Comparing Fig.(|T3|a) and Fig. (|T5|b) , 
we can see that for smaller values of x, where SC are bigger, the Q 2 dependence of the width 
is stronger, which leads to a broadening of the b t profile. This effect reflects the x and Q 2 
dependence of the mean squared radius of the interaction < b 2 >, defined by 


/ b 2 d 2 b t xG(x, Q 2 ,b t ) f nb 2 db 2 xG(x, Q 2 ,b t ) 
f d 2 b t xG(x, Q 2 , b t ) = xG(x,Q 2 ) 


(33) 


In Fig. pfi]. < b 2 > is plotted for different values of x, Q 2 and R 2 . Due to the SC, < b 2 > 
increases as SC increase, i.e., for x —> 0. In HERA kinematic region, even for A 2 = 5GeV~ 2 , 
< b 2 > increases smoothly. 

This result should be compared with the expectations for the soft Pomeron exchange. 
For the soft Pomeron, the slope B is equal to B = B 0 + a' P ln(l/x), where a' P does not 
depend on Q 2 being a' P = 0.25 GeV~ 2 [^]. We parameterize the behaviour of < b 2 > as a 
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Figure 12: Ri( Eq. (fTOJ), uj (Eq. (fn\)) and^(Eq. (f^J ) function calculated from expression 
with b 2 = 0. The results are plotted as functions of y for different values of Q 2 . 
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straightline for 5 < y < 10, < bf > = 4B(y = 5) + 4 a' e ^ ln(l/x). From Fig.[L4|, we can 
conclude that a' e ff shows a Q 2 dependence unlike for the soft Pomeron and its value is rather 
small both for R 2 = 5 GeV~ 2 and R 2 = 10 GeV~ 2 . Indeed, ~ 0.015 GeV~ 2 « A a' P 

for Q 2 = 2.5 GeV 2 and a! e ^ « 0.028 GeV~ 2 « | for Q 2 = 20 Gel/ 2 and i? 2 = 5 GeV~ 2 . 
This straight line approximation is not valid if we go beyond HERA kinematic region. 


xG(b, 2 )/xG(b 2 =0) 



0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 

b , 2 


xG(b, 2 )/xG(b 2 =0) 



Figure 13: The values of xG(x,Q 2 ,b t ) normalized at bf = 0 as functions of b t for different 
values of x, Q 2 and R 2 = 5 GeV~ 2 . 


2.7 The bt dependence of the eikonal Fo structure function. 


Here we will investigate the x and Q 2 dependence of the slope B as in the previous subsection, 
for the F 2 structure function. We define the modified b t dependent F 2 structure function as 


F 2 (x,Q 2 ,b t ) = Ff c (x,Q 2 ,b t ) 


\x,Q 2 ,b t ) + F 2 D 


0, Q 2 ,bt). 


Each term in the above expression obeys the normalization condition 

F 2 (x, Q 2 ) — j d 2 b t F 2 (x, Q 2 , b t ). 


The first term in Eq. 


F. 


sc 


is the bt dependent Mueller formula 
(x,Q\b t ) = f 3 f Q dQ a { 


(34) 


(35) 


(36) 
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Figure 14: The mean radius of the interaction < b 2 > for xG (Eq. 


where [[J 

(1(X, Q' 2 , b t ) = Q' 2 )S(b t ) , (37) 

and S(b t ) is the profile function. The Born term for the SC is written as 


F PORN (x, Q 2 , b t ) = £ a s Z 2 f nQ d(lnQ' 2 )xG DGLAP (x , Q' 2 ) 5(6*) . (38) 

97T i 3 J\nQl 

The b t dependent DGLAP evolved F 2 is given by (see Eq. (|28|) ) 

F? GLAP (x, Q 2 , b t ) = F 2 dglap (x , Q 2 )S(b t ), (39) 

where F PGLAP (x,Q 2 ) is calculated from the GRV distributions for valence, sea, and charm 
components as already discussed in the F 2 section (Eq. (|25D ). All the above expressions for 
F 2 fullhll the normalization (|35|). 

In Fig.[T3] we show the b t dependence of F 2 (x } Q 2 ,b t ) normalized to b t = 0 at different 
values of x, Q 2 and R 2 = 5 GeV~ 2 . At low values of x, the SC lead to a broadening of 
the b t profile, but the effect is less pronounced when compared to the b t profile of the gluon 
distribution. 


We calculate also the mean radius < b 2 > for the F 2 structure function, which is given 
by the expression 


Jd6;6 2 F 2 U,Q 2 ,6 2 ) 

f db? F 2 (x,Q*,bf) 


(40) 


25 























F 2 (b t 2 )/F 2 (b t 2 =0) 


F 2 (b, 2 )/F 2 (b, 2 =0) 




0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 


0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 


Figure 15: The values of F 2 [b 2 ) normalized at b 2 = 0. 

The results are plotted in Fig.|TB]. One can see that the effect of the SC on the value of 
< b 2 t > is smaller in the F 2 case then in the xG case. From Fig. [1^ one can calculate 
the value of ot! e ^ which is equal ot! e ^ ~ 0.015 GeV~ 2 ~ Jo a 'p f° r Q 2 = 2.5 GeV 2 and 
a' e ff ~ 0.017 GeV~ 2 for Q 2 = 20 GeV 2 . Comparing with a' e ^ calculated from the gluon 
distribution, we conclude that a' e ^ for F 2 presents a less pronouced Q 2 dependence. 


2.8 Energy sum rules. 

It is well known that the deep inelastic structure functions should obey the energy sum rules 
which look as follows 

[ dx{xT(x,Q 2 ) + xG(x, Q 2 ) } = 1 , (41) 

Jo 

where xG is the gluon structure function and S(x, Q 2 ) = J2f [Qi( x , Q 2 ) + Q 2 ) ] where 

q ( q ) denotes quark ( antiquark) structure function. 

The SC in our approach violate the energy sum rules, since we heavily used in the Mueller 
formula the leading ln{l/x) approximation of perturbative QCD in which the recoiled energy 
was neglected. However, the main contribution to Eq. ([Il]) comes from the value of x which 
are not small (x ~ 0.5) where we expect that the SC are negligible. Since we took into 
account the full DGLAP evolution equations in the first term of the modified Mueller formula 
we expect that this formula will give us a small correction to Eq. fl4lD . In Fig.[l7| we plot the 
energy sum rules as a function of Q 2 with quark and gluon distributions calculated from the 
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Figure 16: The mean radius of the interaction < b\ > calculated using F 2 (Eq. as a 
weight function ( Eq. M»- 


modified Mueller formula. One can see that the corrections to the energy sum rules is really 
small and can be neglected. The corrections go from zero to 2 % as ln(Q 2 /GeV 2 ) goes from 
the initial value ln(Q 0 ) from 1 to 4. 


2.9 Where is the BFKL Pomeron? 


In this section we address the question why we do not see any manifestation of the BFKL 
Pomeron [Q, since the value of the anomalous dimension ( 7 ) has reached the value First, 
let us summarize what we have learned about the BFKL dynamics during the last two 
decades. 


Lin the BFKL equation all terms of the order of (yf)" - have been taken into account. 
They generate the BFKL anomalous dimension of the form [§],p0 


7 


BFKL, 


D) = 


a s N c 1 , 2c4iV c 4 C(3) 1 


7r u 


+ 


+ E *»(-) 


7T 


CcT 


n=5 


U! 




1 

2 + 


u L —u 


A 


(42) 


where r j BFKL (u) = u>i) = From Eq. (42) one can see that the anomalous dimension does 
not exceed much the value 4 and the second term in Eq. ([I2|) gives rise to the diffusion in 
log of transverse momenta. 

2. The Monte Carlo simulation of the BFKL dynamics as well as the semiclassical solution 
of the evolution equation (see Ref. ED shows that we can see the difference between the 
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Figure 17: The energy sum rule calculated from expression Eq. with the quark and 

gluon distribution given by the modified Mueller formula. 


BFKL and the DGLAP evolutions only for the photon virtualities Q 2 which are close to the 
initial one ( Qq). 


3. The BFKL diffusion leads to the breakdown of the operator product expansion [p| 


in 


the region of small x. We can trust the BFKL equation only in a limited region of x, namely 

x 0 


22 


In— < ca s (Q 0 ) 

X 


2\—3 


(43) 

with sufficiently small numeric coefficient c (see Ref. |22| for details ). The attempts to take 
into account the running cig in the BFKL equation give even more restrictive bound pi 
namely 

In— < ca s (Q 0 

x 

Therefore, we expect that the BFKL evolution can be visible in a limited range of x at 

Q 2 - Ql 

On the other hand the SC break the operator product expansion at any value of Q 2 in 
the wide region of x, namely || 


•\ 2 \ —- 
3 


In— > c'a s (Q 2 ) 2 . (44) 

x 

This inequality can be derived from the Mueller formula and corresponds to /c > 1. 

Our answer to the question formulated in the beginning of this section is as follows. The 
SC turn out to be more important than the BFKL resummation in the whole kinematic 













region at x < 10 2 . Actually, one can conclude this just looking at Fig.l and Fig.2, since 
< 7 > —> | in the kinematic region where k > 1 . 


However, to come to more definite conclusion it is very instructive to notice ( see Fig. 
|7j, that the SC in the eikonal approach reduce considerably the value of < 7 >. Even at 
R 2 = 10 GeV~ 2 , < 7 > becomes smaller than \ for any value of x < 10 ~ 2 at Q 2 > 2.5 GeV 2 . 
The situation becomes more pronounced at b t — 0 ( see Fig®. Indeed, for b t — 0 , < 7 > 
turns out to be smaller than | for any value of Q 2 and x. However, we would like to recall 
that our calculation of the SC in Glauber - Mueller approach depends on the hypothesis that 
we made on the large distance contribution to the Mueller formula ( see Eq. (Q)). In Figs. 
[7| and U we calculated the SC originated from distances r±_ < 1 GeV -1 . To illustrate that 
our conclusion on the BFKL contribution does not depend on large distances contribution 
we calculate the function that gives also the average anomalous dimension when the SC are 
small, 

1 dxG MF , A . 

<7>_ xG GRV dln(Q 2 ) ' ( 5) 


The SC contribute to this function only at small distances (r± ~ ( see Eq. ([]) and 

the discussion in section 4.1) and, therefore, we can calculate them without any uncertainty 
related to the unknown “soft” processes. One can see in Figs. |1^ that this function turns 
out to be smaller then |. This fact supports our conclusion that the BFKL Pomeron cannot 
be seen and it is hidden under large SC. We want also to add to the previous discussion that 
we do not believe that R 2 = 10 GeV~ 2 is a correct radius for a proton as we have discussed 
in the introduction. 


It should be stressed that the SC to the deep inelastic structure functions are the smallest 
ones since in other processes we expect smaller values for the radius. For example, in the 
inclusive production of high transverse jet in DIS ( so called hot spots hunting |24|]) with 
transverse momentum of the jet of the order of Q , R m A and we expect that the SC will 
be more important than the BFKL emission. Therefore, we think that the chance to see the 
BFKL Pomeron is ruled out, following this formalism, even in this specially suited process for 
the BFKL contribution. We have to look for other processes to recover the BFKL dynamics. 


3 First corrections to the Glauber-Mueller Approach. 


In this section we discuss the corrections to the Glauber approach (the Mueller formula 
of Eq. (|T7| ). To understand how big could be the corrections to the Glauber approach we 
calculate the second iteration of the Mueller formula of Eq. As has been discussed, 
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Figure 18: The average anomalous dimension for small SC (Eq. &)■ 
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Eq. (|17|) describes the rescattering of the fastest gluon ( gluon - gluon pair ) during the 
passage through a nucleon. In the second iteration we take into account also the rescattering 
of the next to the fastest gluon. This is a well defined task due to the strong ordering in the 
parton fractions of energy in the parton cascade in leading ln{l/x) approximation of pQCD 
that we are dealing with. Namely 

xb < x n < ... < x\ < 1 , (46) 

where 1 corresponds to the fastest parton in the cascade. 

In the first iteration, we evaluate the SC inserting the gluon distribution xG DGLAP (x : Q 2 ) 
into the modified Mueller’s formula, Eq. ([17:). For the second iteration we substitute the 
correction term of the first iteration 


xG = xG 1 (x, Q 2 )(Eq. (17)) - x G dglap (x, Q 2 ) ; (47) 


in Eq. ( |TTp and add the first iteration to the result. 

The substruction of xG DGLAP in Eq. ([47]) means that, doing the second iteration, we take 
into account only two or more interactions of the next to the fastest gluon with the target 
nucleon. We have to do such a substruction to avoid double counting since the one iteration 
of all gluons with the target has been included in the DGLAP evolution or in other words 
in the Born term of our approach ( see Ref.[11] for details). 

Fig|| shows the second iteraction calculations for our three standard observables: R\, 
< lv > and < 7 >. From the figure one can see that the second iteraction changes crucially 
the values of these parameters. As it was already noted for the nuclear case[ll], a remarkable 
feature is the crucial change of the value of the effective power oo(Q 2 ) for the “Pomeron” 
intercept which tends to zero at HERA kinematic region, making possible the matching 
with “soft” high energy phenomenology, even in the nucleon case. It is also important to 
note that the second iteration makes more pronounced all properties of the ratio R\ and the 
anomalous dimension ( < 7 >). 


However, all these features of the second iteration calculations we have to consider with 
many restrictions. Indeed, the only conclusion which we can derive from our study is the fact 
that the correction to the Glauber (eikonal ) approach turns out to be large in the region of 
small x ( x < 10~ 2 ) and becomes of the order of the Erst iteration for x < 1 CP 3 . As it has 
been already discussed in 0 ED, it occurs because the second iteration gives a correction 
of the order of a^ln^Q 2 / Ql)ln(l / x ), which is close to unity for the DLA of pQCD ( for the 
DGLAP evolution at small x). 


Therefore, the iteration procedure cannot be an effective way to solve that problem. We 
have to develop a different technique to take into account rescatterings of all the partons 
in the parton cascade which will be more efficient than the simple iteration procedure for 
Eq. (El- This is the subject of the next section. 
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4 The general approach. 


4.1 Why equation? 


We would like to suggest a new approach based on the new evolution equation to sum all 
SC. However, first of all we want to argue why an equation is better than any iteration 
procedure. To illustrate this point of view let us differentiate the Mueller ( see Eq. @ ) 
formula with respect to y — ln(l/x) and £ = InQ 2 . It is easy to see that this derivative is 
equal to 


d 2 xG(x, Q 2 ) 
dydi 




Q%) | 


(48) 


A nice property of Eq. (|45| ) is that all quantities in Eq. ( |45|) enter at small distances, therefore 
it is under theoretical control. Of course, we cannot avoid all difficulties just changing the 
solution procedure. Indeed, the nonperturbative effects coming from the large distances 
are still important but they are all hidden in the boundary and initial conditions for the 
equation. Therefore, an equation is a good ( correct ) way to separate what we know ( small 
distance contribution) from that we don’t ( large distance contribution). 


4.2 The generalized evolution equation. 


We suggest the following way to take into account the interaction of all partons in a parton 
cascade with the target. Let us differentiate the ^-integrated Mueller formula of Eq. (|TT]) in 
y = ln(l/x) and £ = 1 ii(Q 2 /Qq). It gives 

a2x g y g^ ] = {C + ln(K G (x',Q 2 )) + Ei(kq(x i , Q 2 ))} . (49) 

where kg{x, Q 2 ) is given by 


DGLAP 


(*,Q 2 ) 


N c a S TT 
2 Q 2 R 2 


xG DGLAP ( x , Q2) 


(50) 


The expression 

9 2 k g (y,£) 

dydl; 


can be rewritten in the form ( for fixed cxs ) 


+ 


dp G (y,0 N c a s 


dy 


7T 


{ c + ln{K G ) + Ei(k g )} = F(kg) ■ 


(51) 


Now, let us consider the expression of Eq. (|5l|) as an equation for k G - This equation has 
the following nice properties: 
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1. It sums all contributions of the order (as2 /£) n absorbing them in xG(y , £), as well 
as all contributions of the order of K n . Therefore, this equation solves the old problem, 
formulated in Ref. ||, and for N c —> cx) Eq. O gives the complete solution to our problem, 
summing all SC; 

2 .The solution of this equation matches with the solution of the DGLAP evolution 
equation in the DLA of perturbative QCD at k —> 0; 

3. At small values of n ( n < 1 ) Eq. ([5l]) gives the GLR equation. Indeed, for small k 
we can expand the r.h.s of Eq. ([E]) keeping only the second term. Rewriting the equation 
through the gluon structure function we have 


d 2 xG(y,£) 

dydS, 


as N c 

7T 


xG(x, Q 2 ) 


ajN 2 c 

8Q 2 R 2 


(xG(x,Q 2 )) 2 , 


(52) 


which is the GLR equation [5|] with the coefficient in front of the second term calculated by 
Mueller and Qiu 0. 

4. For asyC, ~ 1 this equation gives the Glauber ( Mueller ) formula, that we have 
discussed in detail. 


Therefore, the great advantage of this equation in comparison with the GLR one is the 
fact that it describes the region of large k and provides the correct matching both with the 
GLR equation and with the Glauber ( Mueller ) formula. 

We propose also an equation for the ^-dependent gluon distribution. Let us consider the 
^-dependent Mueller formula Eq. (|3l|) . Differentiating this expression in y = In(1 /re) and 
£ = ln(Q 2 /Qg) we obtain 


d 2 xG(y, £, b t ) 
dyd £ 

where Kq(x, Q 2 1 b t ) is given by 


..DGLAP 


(x,Q 2 ,b t ) 



e ~K, G (x,Q 2 ,bt) 


} 




N c a s n 2 
2Q 2 


xG dglap 


(x,Q 2 ,b t ) . 


(53) 


(54) 


As it was discussed before, expression ( |53| ) can be rewritten in the form (for fixed ag) 

\x,Q 2 ,b t ) } = p^DGLAPj _ ( 55 ) 


d 2 KG{y,£,jbt) + dhi G (y, £, b t ) _ N c a s ^ 


d\M 


dy 


_^DGLAP („ r >2 


7T 


Now, let us consider the expression (|55D as the equation for kq. It means that we propose 
the following equation 


d 2 K G (y,£,b t ) + dK G {y,t,b t ) 


N c a s 


dydi 


dy 


71 


{ 1 - } = Fb (K G (x,Q 2 ,b t )) . (56) 


34 















This expression is an equation for n(y,£,b t ). The ^-dependent gluon distribution can be 
obtained from the solution of the above equation using the expression fl54l). 

We claim that Eq. (|5l|) and Eq. ([5(|) solve the problem of summing all Feynman diagrams 
in the following kinematic region: 


a s ln(l/ x) ln(Q 2 /Qo) 


as ln(l/x) < 1 , 


a s hi(Q 2 / Q\ 
as k < 1 . 


< 1 


as -C 1, 


(57) 

(58) 

(59) 


The first two equations specify that we are doing our calculations in the double log approx¬ 
imation of pQCD while the third one introduces the new parameter n (see Ref. |5[) which 
takes into account the parton-parton interaction inside of the perturbative cascade. To il¬ 
lustrate how this set of parameters works let us consider an example simpler than the deep 
inelastic scattering with the proton, namely, the deep inelastic scattering with the meson 

The mass of heavy quarks (?7Tq) sets the scale for the 
oc - 1 7 . Therefore, we can safely apply 

m Q a S( m g) ’ J 


built with heavy quarks (onia) 
distance in this interaction, and we have r± 

pQCD for this problem. 


At order ag(mn) of pQCD, using the set of parameters mentioned above, the DIS can 
be described by the diagrams in Fig.|20|a. The first diagram corresponds to the DGLAP 
evolution equations while the second one is the Glauber - Mueller rescattering in the Born 
approximation. Going to the next order of pQCD, one can see that the emission of one extra 
gluon can be reduced: (i) to the next order correction of the DGLAP evolution equations 
(Fig.|2U|b); (ii) to the emission of the extra gluon with the rapidity close to the photon and 
its interaction with the target ( see Fig.^Olc ) which has been taken into account in the 
generalized Eq. ( 0 ) and (ii) the emission of the extra gluon with the rapidity close to the 
target one ( see Fig.pl]a). The last one is the first from the more general class pictured in 
FigPTIb ( so called enhanced diagrams) that has not been included in Eq. (|5T|). However, in 
Ref. |J it has been shown that integration over k 2 as well as over rapidity y' ( see Fig.|2T|b 
) concentrated at y' close to the target rapidity and they give a small contribution of the 
order as in comparison with the diagrams shown in Fig.pOl. Therefore, we can neglected 


such contribution. However, there is a danger in such an assumption, because it was made 
relying on perturbative estimates while Eq. (|5l|) is a nonperturbative one. Nevertheless, 
the enhanced diagrams reveal themselves a class of diagrams with quite a different topology 
which should be discussed and calculated separately. We are going to do this job in our 
further publications. We must admit that only such study will be able to establish the real 
kinematic region where we can trust Eq. 0). 
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a) b) 




Figure 21: The enhanced diagrams of the deep inelastic scatering with heavy quarks (onia) 
(not included in Eq. (jdlj)). 


Eq. © is the second order differential equation in partial derivatives and we need 
two initial ( boundary ) conditions to specify the solution. The first one is, at fixed y and 

Q 2 —»■ oo, 


2 Q 2 R 2 


The second one we can fix in the following way: at x = x$ (y — yf) which is small, namely, 
in the kinematic region where a^yf < 1 


K 


Hi 


N c a S TT 
2 Q 2 R 2 


xG(x, Q 2 ) , 


(60) 


where xG is given by the modified Mueller formula ( see Eq. (p~7|)). Practically, we can take 
Xq = 10 -2 , since the corrections to the MF are small at this value of x — xq. However, as we 
have discussed, we cannot believe in any treatment of the SC at large distances. Therefore, 
we will take Q g = 1 GeV 2 as the initial value of Q 2 in the modified Mueller formula. Thus, 
the boundary condition (|R|) will be considered only for Q 2 > Ql = 1 GeV 2 . For Q 2 < Ql 
we have to specify the initial conditions and we will discuss this point later. 

The boundary condition for the ^-dependent equation (^6|) is fixed in a similar form, 
namely, we suggest the following initial condition for Q 2 > Qq where the value of h at 
x = x 0 calculated from the b t dependent gluon distribution 
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K 


K i: 


N c a s vr 2 

2 Q 2 


xG(x,Q 2 ,b t ) . 


(61) 


4.3 The asymptotic solution. 


First observation is the fact that Eq. (^l|) has a solution which depends only on y. Indeed, 
one can check that k = K a sym P (y ) is the solution of the following equation 


dn, 


asymp 


dy 


T( K asymp ) 


The solution to the above equation is 

r fcasymp (y) dn' 


o(y=yo) F(k') 


= 2/ — 2/o 


(62) 


(63) 


It is easy to find the behavior of the solution to Eq. ( |63| ) at large value of y since 
F{k) —> dg In k at large k ( dg = — ag ). It gives 




asymp 


At small value of y, F(k) 


k, 


a s y In (a s y) at a s y > 1 
dgK and we have 

(:y = vo) . 


asymp 


ft, 


asymp \ 


(64) 


(65) 


We claim this solution is the asymptotic solution to Eq. (|5iP. To prove this we have 
to consider the stability of the asymptotic solution. It means, that we solve our general 
equation looking for the solution in the form 


ft ft asymp ( V) + AK(y,f-f 0 ) , 


( 66 ) 


where Ak is small (Ak <C n a sym P ) but an arbitrary function at £ = £ 0 - We have to prove 
that Eq. (pip will not lead to big Ak (Ak K asymp ) at large £. 

The following linear equation can be written for Ak 


d 2 A K(y,£) dA k(?/,£) 
dy <9£ dy 


dF{n) 

dn 


R —fcasymp 


(y) Ak(j/,£) . 


(67) 


In Ref. |Tf was proven, that the solution of Eq. (|67|) is much smaller than k, since dF J^ —■> 0 
at large y. 
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The ^-dependent equation also has an asymptotic solution. Indeed, one can check that 
k = K aS ymp{y,bt) is the solution of the following equation 


which can be rewritten as 


dn, 


asymp _ w , 

-T h+ ( fci 


dy 


bt \ ™asymp ) > 


o(y,bt) 


dn' 


?{y=yoM) Fb t (n') 


= y-y o • 


( 68 ) 


(69) 


This equation has the solution 

Kasym P (y,bt) = ln{l + - 1)/^"^}. (70) 

It is easy to find from Eq. ( |70| ) the behavior of the solution of Eq. ( |68| ) at large values of y. 
It gives 

b^asympiy i bt) * 1/ clt 1 • (*71) 

At small values of y, we have 

Kasymp(y,b t ) —> K asymp (y = y 0l b t ) e as(y ~ yo) . (72) 


Therefore, the asymptotic solution has a chance to be the solution of equation (|T]) (fl56|) 
for ^-dependent case) in the region of very small x. To prove that, we need to solve the 
equation in the wide kinematic region starting with the initial condition. We managed to 
do this only in the semiclassical approach. 


4.4 The semiclassical approach. 


In the semiclassical approach the solution of eq. (^Tj) is supposed to be in the form 

k = e s , (73) 


where S' is a function with partial derivatives: ^ = u and f| = 7 which are smooth 
functions of y and £. It means that 


d 2 S dSdS 


d£dy dy 

Using eq. 0). one can rewrite eq.([H]) in the form 

dsos os , s . 
ei di + 8» = e F(e ) s 4(S) 


(74) 


(75) 
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or 


a ;(7 + 1) = $(S) . 


(76) 


For an equation in the form 


F{€,y,s, j,w) = 0 , 


(77) 


we can introduce the set of characteristic lines (C(y), S(y),u(y),'y(y)), which satisfy a set 
of well defined equations (see, for example, Refs. |3l|], [Q for the method). In Ref. |ll 
we developed a detailed calculation of the solution. Here, we will concentrate on the b t 
dependent evolution. To solve the b t dependent equation in the semiclassical approach, we 
need only to substitute the nonlinear term F bt (K), instead of F(k), in the expression (|75|). 
Using eq. (|75D and eq. (j76|) , written for we obtain the following set of equations for the 
characteristics 


d£ _ $ h (S) dS _ 27 + I ( tt = c f>' 1 

dy (7 + 1 ) 2 ’ dy (7 + I ) 2 bt{ h dy SA 1 + l 


(78) 


where $' S;6t = The initial conditions for this set of equations we derive from eq. (Hi. 

being 


So 

7o 


lnK,in{yo, £ 0 ) > 


dlnhii n (yo, £) 




«=So 


(79) 


The DGLAP equation for k is obtained taking the term in the curly in r.li.s of Eq. O 
equal to «, which gives and = 0 in Eq. ([7^). The main properties of these 

equations have been considered in Ref. [IT] analytically. Here, however, we restrict ourselves 
mostly to the numeric solution of these equations. 


For Q 2 < Qq we cannot specify the initial condition. To find characteristics which 
correspond to such small virtualities, we put a boundary condition at Q 2 = Q q, namely 
Abound = k(v,£ = &)• However, this boundary condition cannot be arbitrary at x — > 0, but 
it should satisfy the equation. On the other hand we wish to incorporate in this boundary 
condition everything that has been known on the parton densities experimentally. We suggest 
the following boundary condition: 


Abound — k grv (Q 2 — Ql,x;b t ) at x > x 0 ; 


(80) 


Abound 


asymp 


at x < x 0 , 
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where the value of xq we find from the equation 


k grv (Q 2 = Qq, x = x 0 ; b t ) = 1 . 


Using Eq. (80) we can specify the boundary condition for our set of equations, namely 


*So(£ £cb bt) ^(^boundiyOi £o; bt) , 

/ c c \ dXwiyKbound^QiVibt) , 

— qo) — \y=yo 


( 81 ) 


Using the boundary conditions of Eq. (81), it is better to rewrite Eq. 


in the form: 


dy 


CO 


dS 


®b t (S) ’ dy 


= 2u; + 


u 


div 


®b t {S) ’ dy 


= u 


VsaW 

$ bt (S) 


(82) 


We would like to emphasize that our suggestion ( see Eq. (g0|) ) is still rather arbitrary since 
the equation for xq has no deep theoretical arguments. Actually, what we know theoretically 
is the fact that K boun d —> K aS ymp at k grv > 1. This is the reason why we decided to restrict 
ourselves to the solution of the equation for the trajectories that started at y = 5 with the 
values of Q 2 > Q d = 1 GeV 2 . We are striding other trajectories on a separate paper. 

For the numerical solution we use the 4th order Runge - Kutta method to solve our 
set of equations with the initial distributions of Eq. m- For a short notation we will use 


AGL to our nonlinear evolution equation Eq. (51) and AGLb t for the bt dependent nonlinear 
evolution equation Eq. (^). Fig.|22| shows the general features of the solution. Fig.f^a shows 
the characteristic curves in the y v.s. £ plane for different initial conditions. Fig.|~ 


shows 

the evolution of the respective 7 values with y. The initial conditions, plotted in Fig.(^3|), 
are calculated for y = 4.6 (x ~ 1CT 2 ) and Ql from 1 GeV 2 to 3.4 GeV 2 with Ft 2 = 5 GeV~ 2 . 
When 7 goes to zero as y grows, the nonlinear effects play an important role. The respective 
trajectory tends to a vertical line, and the AGL solution tends to the asymptotic one. When 
7 goes to a constant, the AGL solution tends to the DGLAP one. 


4.5 Comparison between different approaches to the SC 

Using the characteristics and the initial conditions we are able to find the solution of AGL 
equation as a function of y and £. For a given point (y 1 , £1 ), we vary the initial condition 
until we find the characteristic which passes throught the point. From the characteristic we 
can reconstruct the solution. The same procedure can be used for GLR equation. In Fig.|24| 
we compare the k values obtained from the solution of the nonlinear equation ([5l|) (AGL) 
and the GLR equation. The values of k calculated from the modified Mueller Formula ( 
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Figure 22: The trajectories (a) and the 7 values (b) for the solution of the generalized 
evolution equation for nucleon with R 2 = 5 GeV~ 2 . 
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k values 


k values 





y=ln(1/x) 


Figure 24: The n values solutions of the nonlinear equation (AGL) and the GLR equation. 
The values of k calculated from the modified Mueller Formula for the gluon distribution (MF) 
and from the DGLAP gluon distribution is also shown. 
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MOD MF) and from the DGLAP evolution equation (GRV distribution) are also shown. 
The results are presented as functions of y for several values of Q 2 . The initial condition is 
the same for all equations. 

We can see from Fig|24]that the SC predicted by the MOD MF, the AGL equation and by 
GLR equation have the same order of magnitude at HERA kinematic region (5 < y < 10). 
For very small x (y > 10) the results are quite different. The values of n obtained from 
the AGL equation suffer a sizeable reduction when compared with those obtained from the 
modified Mueller formula. The AGL solution increases monotonically with y, but much less 
then the Mueller formula. The GLR equation gives much stronger SC. In fact, the GLR 
equation predicts saturation at very small x values. 

In Fig.^, we present the solution of the AGL equation compared with the AGL asymp 
results. For each value of Q 2 , the asymptotic equation is solved backwards from the value 
of k, solution of AGL equation (we took y = 29.1). We can see from the figure that the 
asymptotic solution crosses the AGL solution twice for Q 2 < 5 GeV 2 . We suppose it occurs 
because the AGL solution have not yet reached its asymptotic value for that y. 


We would like to draw your attention to the fact that our asymptotic solution turns out 
to be quite different from the GLR one. The GLR solution in the region of very small x 
leads to saturation of the gluon density [[TT], |32], pT3| . Saturation means that k tends to a 
constant in the region of small x. The solutions of Eq. ( |5T1) approach the asymptotic solution 
at x —> 0, which does not depend on Q 2 , but exhibits sufficiently strong dependence of k 
on x ( see Fig.|25| ), namely n oc as ln(l/x) lnln(l/x). The absence of saturation does not 
contradict any physics since gluons are bosons and it is possible to have a lot of bosons in 
the same cell of the phase space. We should admit that A. Mueller first came to the same 
conclusion using his formula in Ref. pi 


To complete our discussion, we present in Fig|26] the solutions of the ^-dependent non¬ 
linear equation (AGLfeJ calculated at b 2 = 0 for different values for Q 2 and K 2 = 5 GeV~ 2 . 
One can see that the scale of the SC is much bigger for bt =0 than for integrated over bt 
case. The solutions approach the straight line behaviour for y close to its initial value. It 
occurs because the AGL bt equation reaches its asymptotic behaviour for k( b t ) close to the 
initial condition. Therefore, we can use the asymptotic solution as a good approximation to 
the general solution of the equation at small values of b t . 


4.6 The F 2 structure function from the general approach. 

In order to obtain the F 2 structure function from the general approach, we calculate the 
valence and charm components of F 2 as described in the section [ZT|. We suppose that 
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AGL (b 2 =0) 



y=ln(1/x) 


Figure 26: The hi values solutions of the b t -dependent nonlinear equation (AGL bt ) calculated 
for bj = 0. 
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Figure 27: F 2 calculated from different gluon distributions using the expression for the 
sea component. 


the sea component is generated mainly from the scale violation mechanism. Thus, the sea 
component of F 2 can be calculated by the simplified formula 


iSEA 


Nf 

(x,Q 2 ) = F 2 (x,Ql) + —Y. 2) 

rl—x 


r\n Q 2 

I In Ql 
X 


a s (Q , Nf)d(\nQ 


/ 2\ 


f ^ dz [x 2 + (1 - z) 2 ]{jf-G(-f—,Q' 2 )} 

JO A. Z A. Z 


(83) 


where F 2 (x, Ql) is the value of the sea component of F 2 at Q 2 = Qo- In this approach, we 
include the SC in the gluon distribution in the expression (|83|). In Fig.^7] we show the results 
of F 2 with the sea component calculated from different gluon distributions using expression 
( |83|) for the sea component. The valence, charm, and the initial value of the sea component 
(at Q 2 = Ql = 1 GeV 2 ) are calculated from the GRV distributions. We can see from the 
figure that the SC predicted by the AGL and GLR equations are negligible for y < 7. For 
y > 7, the SC start to become important. This result can be understood if we look at Fig.|24. 
The n values given by the DGLAP, AGL and GLR equations are very close for Q 2 < 10 GeV 2 
and y < 7, i.e., negligible SC are predicted in this region. The results presented in Fig.|27| 
are consistent with that presented in FigJTl|. 


5 Conclusions. 

In this paper we re-analyse the situation of the shadowing corrections in QCD for proton 
deep inelastic structure functions. As have been discussed in the introduction we view the 
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physical interpretation of HERA data as controversial. Here, we are going to summarize the 
result of our analysis of the value of the SC, listing the questions that we asked ourselves 
and the answers we got in the paper. 

1. Can we provide reliable estimates for the value of the SC in QCD? 


Our answer is yes and, to answer this question, we reconsidered the Glauber - Mueller 
approach for the SC in DIS JTO] and suggested the modified Mueller formula ( see Eq. (|T7|) ) 


to estimate the scale of the SC. We consider these estimates as estimates from below for two 
reasons: (i) in the eikonal (Glauber - Mueller ) approach only the fastest partons interact with 
the target and, therefore, the real value of the SC should be bigger; (ii) we cannot trust the 
pQCD approach at large distances which enter into our formulae, in spite of the fact that the 
Mueller formula gives the infrared stable answer for the deep inelastic structure functions. We 
introduce a scale = -f? and calculate only the SC from the distances r± < r 0 . We chose 
Q o = l GeV, relying on the fact, that all available parameterizations of the deep inelastic 
structure functions describes the data down to photon virtuality Q 2 = Q (j = 1 GeV 2 . We 
checked that the modified Mueller formula gives the small corrections to F 2 (x,Q 2 ) at large 
values of Q 2 as well as to the momentum sum rules. We believe that we are able to give a 
reliable estimate of the minimal value of the possible SC to the parton densities. 


2. What is the value of the SC? 


Using the modified Mueller formula we obtained the estimates for the value of the SC 
given in Figs.(^), (]§!) and ©■ These figures show that the SC are essential in HERA 
kinematic region to gluon densities, especially to < 7 > and < u> >. We calculated also the 
SC for b t distributions and showed that the SC at b t = 0 turn out to be much bigger than 
the integrated over the b t parton densities. One can see from Figs.([12|) that < 7 (b t = 0) > 
is smaller that 1 — \ for Q 2 smaller then 10 GeV 2 and < uj{b t = 0) > reaches the value of 
the intercept of the soft Pomeron making possible the matching of soft and hard processes. 
Therefore, our answer to the question formulated in the title of this subsection is that the 
value of the SC is big enough to be seen in the gluon structure function at HERA kinematic 
region. However, the SC to F 2 (x, Q 2 ) turns out to be rather small. This is a reason why the 
SC has not been observed and the whole issue of the SC has been hidden in uncertainties of 
the value of the gluon structure function in current attempts to describe the HERA data in 
the framework of the DGLAP evolution. 


We would like to point out that all numbers were calculated using the GRV param¬ 
eterization for the gluon structure function and our reasoning for the choice of the GRV 
parameterization is (i) it describes the experimental data for Q 2 > 1 GeV 2 ] (ii) the initial 
virtuality in the GRV parametrization is very low ( Qq = 0.3GeV 2 ) and, therefore, the DLA 
of pQCD works better than in the other parameterizations ( we recall that the Mueller for¬ 
mula for the SC was proven in the DLA). However, we are planning to calculate the value 
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of the SC using other parameterization in future. It remains also for further publications 
the important question how the value of the SC depends on the choice of the initial parton 
distribution at fixed value of x. 

3. Where is the BFKL Pomeron? 


Our answer is that the BFKL Pomeron is hidden under the SC. We argue this point in 
subsection 2.9. Here, we want to repeat that this statement depends crucially on the choice 
of the gluon correlation radius R 2 and we argued that HERA data gives R 2 = 5GeV~ 2 . 
This result means that the break of the operator product expansion appears rather as the 
result of the SC than the BFKL dynamics ( see Ref. ). We also argued that the BFKL 
Pomeron should not be seen in all other processes even in the inclusive production of the 
jet in DIS which was specially proposed to test the BFKL Pomeron |24| | . since for all these 
processes R 2 is smaller than for the proton target. We have to search for new processes to 
recover the BFKL dynamics. It is interesting to notice that the fact that the BFKL Pomeron 
is hidden under the SC is seen directly from Figs.[T] and by comparing two curves: k = 1 
and < 7 > = |. One can notice that the Erst curve is to the right of the second one. It 
means that decreasing the value of Q 2 at any fixed x we meet first the curve k = 1 where 
the SC become essential and only for smaller value of Q 2 we approach the curve < 7 > = | 
where the BFKL dynamics could enter in the game. 

4- How well works the Glauber - Mueller approach? 


We showed that the corrections to the Glauber - Mueller approach are big and it can 
give only estimates from below for the value of the SC ( see section 3 ). 

5. Can we develop a general approach in QCD to calculate the SC? 


Our answer is that we can. We suggest the new evolution equation and discussed its 
theoretical accuracy. We argued that this equation solves the problem of the SC summing 
up all diagrams of the order (ctg ln(l/:r) \\\{Q 2 /Q‘l)) n and (k) u considering agln(l/a;) -C 1, 
as \n.(Q 2 /Ql) C 1, as < 1 and k -C 1, where n = xG(x, Q 2 ). This equation allows 
us to calculate the value of the SC in a wider kinematic region than the GLR equation and 
study the parton density to the left of the critical line of the GLR equation ( < 7 >= 1/2). 

6. Do the parton densities reach saturation? 


Our answer is no. We solve the new evolution equation and found out that the gluon 
density xG(x,Q 2 ) tends to the asymptotic behaviour ( xG(x,Q 2 ) —> xG asym (x } Q 2 ) = 
2R 2 Q 2 ln(l/a;) lnln(l/a;)). We discuss the more detailed comparison of the solution to the 
new evolution equation with the GLR one in subsection 4.3. 


7. Can we match the deep inelastic cross section with the cross section of the real photo¬ 
production? 
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We proved that the solution of the new evolution equation gives the mild behaviour of 
the parton densities with energy (asymptotically, they tend to ln(W) ). It gives us a hope 
to describe the matching between soft and hard processes taken into account the SC using 
the new evolution equation. We suppose to do this in our further publications. 


8. What can be a practical use of our approach? 


We were impressed by the fact that HERA data on F 2 (x, Q 2 ) can be fitted by the simple 
formula jf54]|: 

rp f n 2n i Q 2 i xo 

F 2 {x, Q ) = m log —2 log — , 


Q 


X 


where m = 0.377, xq = 0.074 and Ql = 0.5 GeV 2 . Such a parameterization leads to the gluon 
structure function |34|] xG(x,Q 2 ) = 31og(a; 0 /x), which energy (x) dependence is just the 


same as for the solution to our evolution equation. Therefore, we can interpret our results 
as the theoretical justification of the Buchmuller - Haidt parameterization. Our physical 
idea sounds as follows: everything has happened at x of the order of 10 -2 , at such values 
of x the parton densities have reached their asymptotic values and we see only asymptotic 
behaviour of the parton densities at smaller values of x. We suppose to work out this idea 
in the nearest future. 


9. Has everything been done? 

Of course, the answer is no. We have a lot of problems to solve. First, as we have 
mentioned, we need to estimate the values of the SC using other parameterizations of the 
parton densities but not only the GRV one. We have to study more carefully the dependence 
of the value of the SC on the initial parton distributions and to check how stable are our 
results. Second, we have to find the complete solution for the running ag. The experience 
with the GLR equation shows that the running ag leads to important properties of the SC 
as the appearance of the critical line and so on. Third, we need to take into account that 
the partons in the parton cascade can interact not only with the target but with the other 
partons (so called enhanced diagrams). We can End out the region of the applicability of 
the new equation only after solving this problem. We suppose to do this as the next step in 
our approach. 

We believe that our paper is only the first step in reconsidering the whole issue of the SC 
in DIS. Actually, we showed here, that the new parameter that governs the size of the SC 
is k and the value of the SC is not negligible in HERA kinematic region. This conclusion 
makes rather suspicious the current interpretation of HERA data on the basis of the DGLAP 
evolution equations without the SC and much work is needed to clarify physics of DIS that 
is behind HERA data. 


50 





While completing this paper, our attention was drawn to three recent papers |I6| |37 
which deal with the Glauber - Mueller approach to the SC for the deep inelastic structure 
function Fo(x,Q 2 ) for nuclei and for the hard diffractive production. The results obtained 
in these paper are close to ours. 
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